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DYNAMICS OF A STOCHASTIC RATIO-DEPENDENT 
PREDATOR-PREY MODEL 

NGUYEN THI HOAI LINK, TA VIET TON* 


Abstract. In this paper, we consider a stochastic ratio-dependent predator- 
prey model. We firstly prove the existence, uniqueness and positivity of the 
solutions. Then, the boundedness of moments of population are studied. Fi¬ 
nally, we show the upper-growth rates and exponential death rates of popula¬ 
tion under some conditions. 


1. Introduction 


The stable predator-prey Lotka-Volterra equation 

\xt = (ai - ciyt - bixt)xt, 

\yt = (-02 -I- C2Xt - b2yt)yt, 


( 1 . 1 ) 


has attracted much interest. In this model, x* and yt represent the population 
density of prey species and predator species at time t, respectively; oi is the intrinsic 
growth rate of prey in the absence of predator, 02 is the death rate of predator in 
the absence of prey; bi measures the inhibiting effect of environment on two species; 
Ci{i = 1,2) are coefficients of the effect of a species on the other and those coefficients 
are positive constants. It is well-known that the solution of (HID is asymptotically 
stable. 

For the stochastic predator-prey Lotka-Volterra equation, we have to mention 
one of the first attempts in this direction, a very interesting paper of Arnold et al. 
[2] where the authors used the theory of Brownian motion and the related white 
noise models to study the sample paths of the equation 


dxt = (oi - Ciyt - bixt)xtdt + axtdwt, 
dyt = (-02 -f C2Xt - b2yt)ytdt + pytdwt, 


( 1 . 2 ) 


where {wt,t > 0} is a 1 - dimensional Brownian motion defined on a complete 
probability space with filtration (11, {J^t}t>o, P) satisfying the usual conditions 

(see 0). The positive numbers cr, p, respectively, are the coefficients of the effect of 
environmental stochastic perturbation on the prey and on the predator population. 
In this model, the random factor makes influences on the intrinsic growth rates of 
prey and predator. 

For both models (HU) and H2D, the predator consumes the prey with functional 
response of type CiXtyt and contributes to its growth rate C2Xtyt- For such models. 
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we refer to lai. However, some biologists have argued that in many situations, 
especially when predators have to search for food, the functional response should 
depend on both prey’s and predator’s densities; for example, a Holling’s type II 
functional response in a type of [9l [10] , a ratio-dependent one in a type of 

xt+-yyt 13' model in those papers are deterministic. 

In this paper, we consider a ratio-dependent predator-prey model in the random 
environment in which the white noise makes the effect on both the growth rates 
of species and the inhibiting effects of environment on two species, i.e., oi —oi -I- 
(TiWt, —02 —02 + piWt, —hi —>■ —hi + (JiWt, —1)2 —1)2 + P 2 Wt- The model is as 

follows 


dxt 

< 

dyt 


ai(t) - bi{t)xt - ^tdt 

Xt + e{t)yt_ 

- [cri(t) -I- a2{t)xt]xtdwt, 

-02(t) - b2{t)yt + 

xt + e{t)yt 

- [pi{t) + P2{t)yt]ytdwt. 


ytdt 


(1.3) 


The white noise here is effected by one-dimensional Brownian motion {wt,t > 0} 
dehned on the complete probability space (H, Jt, P) with filtration satis¬ 

fying the usual conditions. Throughout this paper, we suppose that all coefficients 
of deterministic part of m are bounded by some positive constants, i.e., their 
inhmum and supremum on f > 0 are positive, and assume that at, pi{i = 1,2) are 
nonnegative bounded functions. 

The hrst aim of this paper is to study the existence, uniqueness and positivity 
of solutions of system m by using some Lyapunov functions. The next aim is 
investigating the boundedness of moments. And the final one is to show some 
upper-growth rates as well as showing the exponential death rates of species under 
some conditions. 

The paper is organized as follows. In Section 2, we prove the existence, unique¬ 
ness and positivity of solutions. Section 3 produces asymptotic estimations of mo¬ 
ments. In section 4, the upper-growth rates and exponential death rates of species 
are studied. 


2. Global solution 

In this section, we show that the solution of system m is positive and global. 
Denote M" = {(xi, • • • ,a;„) G > 0 (z > 1)}, ]R"o = {{xi, - ■ ■ ,Xn) G 

]R":a;i > 0 (z > 1)} (rz = 1,2). Let g{t) be a function, for a brevity, instead of 
writing g{t) we write g. If g is bounded continuous function on K+o, we denote 

= sup g(t), g^ = inf g{t). 

zgR+o teR+o 

We have the following theorem. 

Theorem 2.1. For any given initial value {xo,yo) G there is a unique solution 
{xt,yt) to (1131) for t > 0. Further, with probability one, is positively invariant 
for (11.311 . i.e., {xt,yt) G a.s. for all t > 0, if{xo,yo) G M+. 

Proof. Since some coefficients of (HH) are not locally Lipschitz continuous, we can 
not say that there is a unique local solution (xt,yt) to (11.311 . However, we consider 












A STOCHASTIC RATIO-DEPENDENT PREDATOR-PREY MODEL 


3 


the following system 


d^t 


< 

drjt 


ai - -(cTi + (J2 exp{^t})^ - hi exp{^t} 
+ (cti -I- (72 exp{^t})dwt, 


Cl exp{?7t} 
exp{6} + e expiry*} 


dt 


-02 - 2 (Pi + P2 expiry*})^ - 62 expiry*} 
+ (Pi + P2 exp|ry*})drt;*, 


C 2 exp|g*} 
exp|6} + e expiry*} 


dt 


( 2 . 1 ) 

with an initial value (^O) Vo) = (lna;o, Inr/o)- Since the coefficients of (12.1|) are locally 
Lipschitz continuous, there is a unique local solution (^*,ry*) to (12.11) for t G [0,t), 
where r is the explosion time (Arnold [1] or Friedman i)- Therefore, by Ito’s 
formula, (xt,yt) = (exp|^*},exp|ry*}) is the unique positive local solution to (11.31) 
for t G [0, t) with the initial value (aio, j/o)- To show the solution is global, we need 
to show that r = 00 a.s. We use the technique of localization dealt with in [3 E]. 
Let fco > 0 be sufficiently large for xq and yo lying within the interval [^, fco]. Let 
us define a sequence of stopping times [S] Problem 2.7, p.7] for each integer k > ko 
by 


Tfc = inf 


t>0:xt^ {y,k) or ry* ^ (ifc) 
k k 


(with the convention inf0 = 00 ). Since Tk is nondecreasing as k ^ cjo, there exists 
the limit Too = limfc_*ooTfe. Then Too < t a.s. Now, we will show that Too = 00 
a.s. If this statement is false, then there exist T > 0 and e G (0,1) such that 
1 P|too <T} > e. Thus, by denoting flk = {xk <T}, there exists fci > /cq such that 


P(r2*j) > £ for all k> ki. 


( 2 . 2 ) 


Let 0i G (0,1) {i = 1, 2). We consider the following function 

2 

V{x,y) = x^^ — In a: -I- — In 1 / — ^ ^ ' 

i=i 

Because x^^ — In a; — > 0 for all x > 0, we have V G (^^(K^, R+o). If 

{xt,yt) G K+, by using Ito’s formula, we get 

dV{xt,yt) = f{xt,yt,t)dt + g{xt,yt,t)dwt, (2.3) 


where 

9{X, y, t) =(01X^1 - l)(tTi + a2x) + (02^’' - l)(pi + P2y), 


fix,y,t) =(01X^1 


1 ) 


oi — bix 


ciy 

X + ey 


+ (^ 2 /^ - 

+ 21^2(02 


1 ) 


-02 - b2y - 


C 2 X 


X + ey_ 
l)x®^ + l]((Ti + 0 - 22 ;)^ 

l)y®" + l](pi +p2y)^. 


(2.4) 


It is easy to see from 0* G (0,1) and from (12.4p that the function /(x, y, t) is bounded 
above, say by M, in x R+q. It then follows from {xt/\Tk, yt/\Tk) G and from 
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(j2.3[l that 

rTATk pTATk /-TATfc 

/ dV{xt,yt)< / Mdt+ / g{xt,yt,t)dwt. 

Jo Jo Jo 

Taking expectations yields 

EV{xTATk,yTArJ < V{xo,yo) + ME{T Ar^) < V{xo,yo) + MT. (2.5) 

On the other hand, for every lu € flk, either Xr^i^) or yr^i^) belongs to the set 
{k, -1}. Then 

1 + In 


y{xTATki^)^yTArki^)) >min -\nk- 


e. 


1 1 1 + In 

w.-'-’-k—— 

= miii(t*‘ - Ink - ^ + 


e,. 


1 1 + In 

+ »— (‘ = ‘'2 


We therefore get from (12.21) that 

EV{xTATk, yrATk ) V(xtat^ , Vtatu )] 


lin I 


>e min < — In /c — 


1 + In 


0,: 


1 1 + In 0* 

_ + i„t — 


(i = l,2)}. 


It then follows from (j2.5ll that 
V (xo jyo) + MT > £ min | fc®* — In /c — 


1 + In 9i 1 


9 . + 9 , 

Letting k ^ oo leads to oo > V{xo,yo) + MT = oo. This is a contradiction. 
Therefore Too = oo a.s. Then r = oo a.s., and {xt,yt) S a.s. The proof is 

complete. □ 


1 + In 


(* = 1 , 2 )}. 


3. Boundedness of moments 

In this section, we shall show the boundedness of species’s quantity moments. 
In order to capture them, we shall distinguish between the following hypotheses of 
random perturbation. 

(HI) The random factor make only effects on the growth rate of population, i.e., 
o’ 2(0 = P 2 {t) = 0 for alH > 0 and > 0; 

(H2) The random factor make effects on both the growth rate of population and 
the inhibiting effects of environment, i.e., (j\, p\ > 0 (i = 1,2). 

Consider 

1 1 

LV{x,y) =-[a,{t) + a 2 {t)xfx^^ + -[p,[t) T P 2 {t)yf y^ 

d'^V 

+ [ai{t) + cr 2 {t)x][pi{t) + p 2 {t)y]xy-^^ 













A STOCHASTIC RATIO-DEPENDENT PREDATOR-PREY MODEL 


5 


the infinitesimal operator of (1131), defined on the space where 

}i(x,y,t) = aix -^- bix , 

a; -I- ey 

J2{x,y,t) = -a2y-\ -^-022/ , 

a; -I- ey 

and for any positive numbers di, 02 , we put 
d2 = min{0,6' {i = 1, 2)}, 0 = ^ , 

+ U2 

di = sup IitTi(t )01 (01 - 1) + ^p?(t) 02(02 - 1) + cri(t)pi(t) 0 i 02 +0iai(t) 

-I-[C2(t) - a2(t)]02|, (3.1) 

Ai = -- (01 + 02){1 — ln(0i + 02 )}, 

M2 

A 2 (a:o, 2 /o) = (0i Inxo + 02 ln 2 /o) + (0i + 02){1 - ln(0i + 02 )} - 

d2 

We have the following theorems. 

Theorem 3.1. Under condition (HI), for any positive 0i,02 and initial value 
( 2 ^ 0 , 210 ) S solution of (11.31) satisfies 

E(a;(^2/f^) ^ exp{Ai -I- A 2 ( 2 ;o, 2 / 0 ) exp{— 1 / 2 ^^}} for all t>0. 

Consequently, limsupj^o,^E(a:(^2/t^) ^ exp{Ai}. 

Proof. Firstly, we prove that 

lK{xf^yf^) < 00 for alH > 0 and 0^ > 0. (3.2) 

Define a function V G C^(]R^,]R+) by V(x,y) = x^^y^^. For any t > 0, using Ito’s 
formula gives that 

dV{xt,yt) = LV{xt,yt)dt + {Oiai + 92Pi)V{xt,yt)dwt. (3.3) 


It is easy to see that 

LV{x,y) = ^ 01(01 - l)cr} + ^ 02(02 - l)p} + 01020 - 1/91 + 01 (mi - 61 a;) 


- 02(02 + 02 //) 


92C2X - 9 iCiy 


Vix,y) 


(3.4) 


X + ey 
<[di - d2{x + y)]V{x,y). 

For every integer k > 1, we define a stopping time = infjt > Q-.xt+yt > k}. 
Then, the sequence {Tk,k > 1} is nondecreasing and by the positive invariance of 
{xt,yt) on R^, we have limfc_>oo Tk = 00 a.s. It then follows from (13. 3 1) that 

l-tATk 

y{xtATk,ytArk) =V{xo,yo) + / LV{xs,ys)ds 

Jo 


rtATk 


[0icri(s) -I- 92Piis)]V{xs,ys)dws- 


Taking expectations of both sides and using (lO) . we have 

rtATk 

EV{xtATk^ytATk) <y{^o,yo)-\-diE / V{xs,ys)ds 

Jo 








6 


NGUYEN THI HOAI LINK, TA VIET TON 


<V{xo,yo)-\-di / EV{XsAT,,,ysArk)ds. 

Jo 


Thus, by using Gronwall’s inequality, we obtain 

EV{xtA.rk,yt/\Tk) < V{xo,yo)exp{dit}. 

Letting fc —)• oo in the latter inequality, we yield that, for all t > 0, lKV{xt,yt) < 
2 /o) expjdit}, from which we deduce (13.21) . 

Next, since V{x,y) = x^^y^'^ < (a:: + we have x -\- y > V^{x,y). It then 

follows from (13.41) that 

LV{x,y) < [di - d 2 V'^{x,y)]V{x,y). (3.5) 

Applying (13.2|) to (0i(l + 0), 02(1 + ^)), we have 

E[V^+^ixt,yt)] =E 

Then, by using Holder’s inequality, yields 

[Ey(xi,y0]'^" <E[V+^(xt,yt)] ■ 

It then follows from (j3.3p and from (13.5p that, for any t > 0 and h > 0, 
EV{xt+h,yt+h) -'^V{xt,yt) < J [diEV{xs,ys) - d 2 EV^^^{xs,ys)] ds 

pt-\-h p 


ei{i+e) 62 ( 1 + 0 ) 
a-t i/i 


< oo 


for all t > 0. 


< J ^diEV{xs,ys) - d2[EV{xs,ys)]^''~^^ds. 

(3.6) 

Putting v{t) = EV{xt,yt), then 0 < v{t) < oo for all t > 0. Further, the continuity 
of v{t) in t can be seen by the continuity of the solution (xt,yt) and the dominated 
convergence theorem. We define the right upper derivative of v(t) by 

From (13.61) . we have 

v{t + h) — v{t) 


- hj. 


h^O 


[(ii'(;(s) — d2n^''"®(s)] ds. 


Letting h —^ 0 gives D^v{t) < v{t)[di — d 2 V^{t)] for all t > 0. Therefore, 


ZI'''[exp{fi2i} Int’(t)] = d2exp{d2t}\nv{t) + exp{d2i} 


D^v{t) 

v{t) 


< 02 exp{02i} Inu(t) +exp{02i}[di — d2V^ {t)] 
= di exp{02i} + 02 exp{02i}[lnt'(t) — v^{t)]. 
It is easy to see that \nx — x^ < —1(1 + ln0) for all a; > 0. Then 


I?^[exp{02t} Inu(t)] < 


01 - -^(1 + In 0)02 

U 


Taking integrations of both sides yields 
exp{02t} Inw(t) <lni;(0) + 

=A2 + Ai exp{02t}. 




exp{02t}. 

[exp{02t} - 1] 
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Consequently, we have lnz;(t) < Ai -I- A 2 exp{—d 2 t}, from which follows the first 
statement of theorem. Letting t —)• 00 in the latter inequality, we get the second 
one. The proof is complete. □ 


Theorem 3.2. Under condition (H2), for any 9i G (0,1], pi G [0)3) and q G 
K+, there exist positive constants Ki = Ki{0i,Q) and K 2 = K 2 [Qi,<ii) {i = 1,2) 
satisfying the following for any initial value (xq, yo) G 


(i) limsupj 

(ii) limsupj 


.E 


Si 


O'. 


<ifi; 


-^00 t 


1 io + <^ 22 /^]ds < d^ 2 - 


Proof. Consider a function Vi R^ defined by Vi(a;, y) = For 

any t > 0, by using Ito’s formula, we have 


dVi{xt,yt) =LVi{xt,yt)dt + |6»i<ji[(Ti(t) + a 2 {t)xt]x^^ 
+ 6*2^2 [pi (0 + P2{t)yt]yt^'^dwt, 


where 


2„6»i 


LVi{x,y) =-9i{0i - l)ci[cri(t) +a 2 (t)xYx 

+ ^02(92 - l)c2[pi(t) + P2{t)yYy^^ 




ai(t) — bi{t)x — 


ci{t)y 


X + e{t)yi 


+ 6*2^2/^ - a2(i) +- b 2 (t)y 


(3.7) 


(3.8) 


X + e{t)y 

Then, from 9 i G (0, 1] and from q G M-i- {i = 1,2), there exists Ki = iFi(0i, 02,‘?i,‘?2) 
such that LVi{x,y) + Vi{x,y) < Ki for all {x,y,t) G R+ x R+q. Applying Ito’s 
formula yields 


d[e‘Vi(a;t,yt)] =e*[yi(a;t,yt) + LVi{xt,yt)]dt 

+ e‘[6»ici(cri + a2Xt)xY + 02<?2(pi + P2yt)yt^)]dwt 
<Kie*dt + e*[6<iCi((Ti + a2Xt)xY + 6'2^2(pi + P2yt)yt^)]dwt. (3.9) 


Using the sequence of stopping times { 7 -^}^^ defined in the proof of Theorem 12.II 
and from (13.91) . we have 

E [e*^'"’’Vi{xtATk,ytATk)] < Viixo,yo) + Ki{Ee*^'">’ - 1) for all t > 0. 

Letting fc —>■ 00 in the latter inequality with a fact that V(xtf^rk i VthTk ) > 0 and 
0 < < e* O.S., and using Fatou’s lemma, we obtain 

e‘EUi(a;t,yt) < Vi{xo,yo) + iFi(e* - 1). 

Therefore, limsup(_^oo EVi(xt, y*) < Ki. 

To prove Part (ii), we consider a function V 2 {x,y) = gix^^ + C 2 y®^- Since pi G 
[0,3), there exist 9i G (0,1) {i = 1,2) such that 0 < pi < 2 + 9i. Then, from (13.81) 
there exists K 2 = K 2 {pi,(,i) such that LVi{x,y) + V 2 {x,y) < K 2 for all {x,y,t) G 
R^ X R-i-o. Using dsii) gives 

Vi{xt,yt) <Vi{xo,yo)-\- / [K 2 -V2ixs,ys)]ds 
Jo 
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+ [ [0l<;i{(7i-\-a2Xs)xl^ -\-02^2{pl-\-P2ys)yt^)]dWs. 

Jo 

Taking expectations of both sides, we obtain 


^Vi{xt,yt) 


^V2ixs,ys)ds < Vi{xo,yo) + K2t, 


from which follows El/ 2 ( 2 ;s, 2 /s)ds < Vi{xo,yo) + K 2 t. Therefore, 


1 /■* 

limsup - / + <r 2 ?/f^]ds < K 2 . 

t->00 t Jq 


4. Upper growth rate estimation 


□ 


In this section, we shall show the upper-growth rates of population under the 
case of the random factor making the effect only on the growth rate of population. 

Theorem 4.1. Under condition (HI), for any 6i>Q and any initial value {xq, yo) S 


In 


lim sup ■ 

t—^oo 


&2 

xt yt 


Int 


— < 6i + 02 


Furthermore, if 9i S [0,1) then for any ci > 0 (i = 1, 2) there exists K = K{9i, q) 
such that 


1 I 

limsup - / (ciccfi -I- C 2 yf^)ds < K 
t—^oo i Jq 


a.s. 


Proof. Firstly, we prove the first inequality. Putting xt = expj^t}, yt = expli^t}, di = 
2 2 ____ 
oi — ^,'d 2 = «2 + ^ and substituting this transformation into system (IE3 , we 



•di - bi exp{^t} - 


ciexp{7?t} 


--02 - h expjiyt} -f 


exp{6} + eexp{? 7 t} 
C2 exp{^t} 


exp{6} + eexp{77t} 


dt + aidwt, 
dt + pidwt, 


(4.1) 


dft = 
dr]t = 


di - hixt - 


ciVt 


xt+eyt 


dt -\- uidwt, 
dt + pidwt- 


[-^2 - b2yt + 

Fix p > 0. Applying Ito’s formula to exp{pt}^t and exp{pt}r]t, from (14.IL we have 


exp{pt}6 =Co 


exp{ps} 


di - biXs - 


ciys 


pexjp{ps}fsds 


/o 


/o 


exp{pt}? 7 t =r]o + / expjps} 
Jo 

rt 


-d2 - b2ys + 

t 


Xs + eys 

(Ti exp{ps}(iws, 

C2Xs 


ds 


Xs -I- eys_ 

+ / pexp{ps}?7sds-I- / pi explpsldtcs- 

Jo Jo 


ds 


(4.2) 


(4.3) 
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We set 


Mit= / aiexp{ps}dws,M 2 t = / pi exp{ps}(iws, 
Jo Jo 


then Mit {i = 1,2) are real valued continuous martingales vanishing at t = 0 with 
quadratic forms 

<Mi,Mi>t= f alexp{2ps}ds, <M 2 ,M 2 >t= f plexp{2ps}ds. 

Jo Jo 

Let £ S (0,1) and 9 > 1. Using the exponential martingale inequality [6l Theorem 
1.7.4], for every k > 1 and i = 1, 2, we have 


P < sup Mjt - - exp{-pfc} < Mi, >t 
Lo<t<fc L z 


> 


6 exp{pfc} 




It then follows from Borel-Cantelli lemma that there exists an C U with P(r2i) = 
1 having the following property. For any w S there exists ki = ki{u}) such that, 
for all k > ki and t G [0, k], 

Mu < ^ exp{-pk} < Ml, Ml >t 

2 £ 

£ r 7 1 /"* 2 rr, 1 7 ^ exp{p/c} 

= - exp{—pk} / ai exp{2ps\ds H- mk, 

2 7o e 

M 2 t < ^ exp{-pfc} < M 2 , M 2 >t infc 

2 £ 


£ / 

= - exp{—pfc} / p\ exp{2ps}(is + 

2 Jo 


9exp{pk} 


In k. 


We therefore have from (14.2p and (1431) that for any w S (~l U 2 and t G [0,k],k > 
where ^ 0 ( 02 ) = ki{ui) A k 2 {ui), 

exp{pt }6 <6 + / pexp{ps}^sds 


exp{ps} 


di - biXs - 


ClVs 


Xs + eps 


ds 


^ /■* 

+ - exp{—pfc} / exp{2ps}(is + 

2 Jo 


9 expjpfc} 


In k 


- f* r , 9exp{pk} /■* 

='S 0 +P / exp|ps}4sas H-lnfc+ / exp{ps| 


ai 


- biXs + ^exp{-p{k - s)}- 

2 Xs + eps-i 


ClPs 


ds. 


(4.4) 


exp{pt}77t <770 + / pexp{ps}risds, 

Jo 

rt 


expjps} 


-■^2 - b2Ps - 


C2Xs 


+ — exp{—pfc} / pf exp{2ps}ds + 

2 Jo 


ds 

Xs + eps_ 

9 exp{pk} 


In k 
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r r 1 J ^ exp{pk} /■* 

=Vo+P / exp{ps}risds -\ -lnfc+ / exp|ps| 


-d2 


- hVs + ^ exp{-p(fc - s)} H- 

2 Xs ^ys - 


ds. 


(4.5) 


From (HU) and (I4.5|) . for any w S n ri 2 and t £ [0, k],k > ko{uj), we have 
exp{pf}( 6 »i^t + 62 ^ 1 ) <( 01^0 + 6 * 2 f?o) + ^ 2 )exp{pfc} 


+ / exp{ps} p( 6 »i^s + 6 > 2 ??s) + 6 ' 1 'di - 6 » 2 d 2 
JO 

e [01(7? + 02P?] r M 

+ —^^ exp{-p(fc - s)} 


- &i0ia:s - &202ys + 


O 2 C 2 XS - 01 Cl 2/s 


ds. 


(4.6) 


Xs + eys 

Since 0 i ,02 G K+, there exists H = H{p, 61 , 92 ) > 0 such that for any {x,y,t) £ 

X M+O; 

92C2X - diciy- 


p( 0 i In a; + 02 In?/) + 0 i 0 i — 02^?2 — 0 i 0 ia: — 02022 / + 


< H. 


X + ey 

It then follows from (|4.6I) that for any w £ fli fl fl 2 and t £ [0, fc], fc > fco(w), 
exp{pt}( 0 i^t + 02 ?/*) <( 0 i 6 + 02 ??o) + expjpfc} 


exp{ps} 


H + 


£{9ia\ + 02 pO 


ds 




1 

V L 


Ff + 


e( 0 iar + 02 /or) 


e 

,u 2 \ 


(exp{pt} - 1 ). 


Thus, 


0 i 6 + ^ 2 ??t <( 0 iCo + 6 ' 2 ??o) exp{-pt} + 


0(01 + 02 ) exp{p(fc - t)} 


In k 


H + 


eiOiaf + 02 pr) 


(1 - exp{-pt}) 


<(«.& + e,m} + '’«’■ + i„ k 


H + 


eiOiaf + 02pf) 


For any w £ fli n fl 2 , k > ko{uj) and t £ [k — l,k], from (14.71) we have 


0 i 6 + 02 T?t ^ 1 


Int 


(9,5.+ + t(h±«£ 2 EM,„t 


ln(fc — 1 ) - 

lf„ e(0iaf + 02pf) 


+ 


P 


{iFT 


(4.7) 


2 
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from which implies 


-f 6 » 2 ? 7 t 6 >( 6 »i -f 6 * 2 ) exp{p} 

lim sup-^- <- 

t^oo mi £ 


Letting e— >1 ,0— and noting P(ni n ^ 2 ) = 1 yields 


lin,snp^l%±^<0i+02 

t^Go ini 


i.e., 


In 


lim sup ■ 

t—^oo 


0 \ 02 
Xt Vt 


Int 


< 0 i + 62 


a.s., 


a.s. 


Now, we prove the remain inequality. Putting V{x,y) = ln(<jia;®"^ +^ 22 /^^) and using 
Ito’s formula, we get easily that 


dV{xt,yt) = 


where 


Putting 


and 


+ <J 22 /^' 

<>20 2 y^^ 


-(oi — bix — 


ciy , 

x + ey^ 


+< 222 /®- 


■(-02 - 622 / H- 


X + ey 

[M0, - l)x<^<{<;,x^< + ^ 22 /^^) - 
2 (^ix^i + ^ 22 /®^)^ 

[<22^2(6*2 - l)2/®"(<2ia:®^ + <222/^^) - <2!6*12/^®^]pi 
2 (cix®i + < 222 /®^)^ 

< 2 i< 226 'i 6 » 2 cripia;®i 2 /®^i 


2 (< 2 ia;®i + ^ 22 /® 0 ^ J 

<2i6»i(Tia;®i + < 226 » 2 Pi 2 /® 


dt 


Cix' 


®i + ^ 22 /®' 


■dwt 


r,r , ft^'ic^ia;®'+ ^ 202 ^ 12 /®" j 

=P{x, y, t)dt H- 5 —- 5 ^^- dwt, 


Pix,y,t) =— 


< 2 iX"i -h < 222 / 

<> 202 y^' 


<>ix^< + <, 2 y^ 


(oi - bix - 


X + ey 


+ 


< 2 ia;®i + < 222 /®- 


■(-02 - 622 / + 


C 2 a: 


+ 


K = 


x + ey 

Q6>icri [<22(6*1 - 1)2/®^ - Qa:®^]a:®^ 
2 (^ix^i + ^ 22 /®^)^ 

^26*2^1 [+(6*2 - l)<r®^ - <222/®^]2/®" 
2 (< 2 ix^i + ^ 22 /®^)^ 
Ci< 226 *i 6 l 2 gipia:®i 2 /®" 

2 (< 2 ix®i + < 222 /®^)^ 


sup [P(a;,2/,6) + (^la:®^ + ^ 22 /®"). 


{x,y,t)GR‘^ xR+o 


Mt = 


<2i6>icria:®i + ^ 26 * 2 Pi 2 /®^ 
< 2 ia;®i + < 22 //®=' 




(4.8) 
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then , t > 0} is a martingale, and by 9i € [0,1), we see that K < oo. The 

quadratic variation of Mt can be shown by using [ 6 j Theorem 5.14, p.25] as follows. 


<M,M >t= 
It is easy to see that 






ds. 


limsup ^ < [max{ 0 i(T“, 6 > 2 CT 2 }]^ . 

t—>oo t 

So, using the strong law of large numbers for martingale [5J Theorem 1.3.4], we 
have 

Mt 


lim — = 0 

t—¥oo t 


a.s. 


(4.9) 


On the other hand, from (14.8p . we have 

Q <V{xt,yt) < [ [if-+<r 2 j/f^)]ds + Mt, 
^0 


from which follows 


j + ‘^'^ys'')ds < if + 


Therefore, by (gH), 


1 i* 

limsup - / (cixfi + C22/f^)ds < if 
t—)-00 t Jq 


a.s. 


□ 


Remark 4.2. For the deterministic version of model dm, i.e., Ui = Pi = Q and 
the other coefficients are constants, it is easy to see that limt_,.oo yt = Q holds under 
some special conditions, i.e., the predator dies out, but it never get limt^oo a;^ = 0 
(if liuit-^oo yt = 0 then liminf^^oo a;* > ^ > 0). However, in the above theorem, 
if K = 0 then both prey and predator die out. This means that a relatively large 
stochastic perturbation can cause the extinction of the population. Further, the prey 
population dies out even if there is no predator and the death rate is so rapid (at 
an exponential rate). We can see that in two following theorems. 


Theorem 4.3. Under condition (HI), if the prey is absent, i.e., Xt = 0 a.s. for all 
t > 0, then the predator dies with probability one. Furthermore, the death rate of 
predator is exponential, i.e.. 


limsup — < - inf 

t—foo i 


02 (t) + 


plity 

2 


a.s. 


Proof. The quantity yt = expjTyt} of predator at the time t satisfies the following 
equation 


dpt = 


Pi 


-02 - y - ^2 exp{77t} 


dt + pidwt. 


Thus, 


Vt =Vo + 


<Po- inf 

i>0 


Pi 


-02 “ y “ &2 expires} 


ds F Mf 


02 (t) + 


plit) 


t + Mt, 


2 


(4.10) 
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where Mt = fg pi{s)dws is a martingale. The quadratic variation of Mt 


<M,M >t= 


[ Plis)ds 
Jo 


satisfying 


]• ^ M,M >t ^2 

lim sup- < Pi . 

1—>oo t 


Using the strong law of large numbers for martingales gives limt^oo -p = 0 a.s. It 
then follows from (14.101) that 


lim sup < - inf 
t —^oo i 




plit) 


and limt^oo yt = 0 a.s. 


□ 


Theorem 4.4. Under condition (HI), if the predator is absent, i.e., yt = 0 a.s. 
for all t > 0, then the quantity of prey satisfies the following 

< 0 then limt_,,oo Xt = 0 a.s. and the prey dies out 

at an exponential rate; 

(ii) //supt>o |ai(t) - = 0 then limt^oo Ex* = 0; 

(hi) 


(i) //supt>o |ai(t) - ^ 


Inxt 

lim sup -- 

t—^oo In t 


< 1 


a.s. 


Proof. Similarly to Theorem l4.31 the quantity Xt = 
the following equation 


d^t 



hexplft} 


expj^t} of prey at time t satisfies 


dt + (Tidwt. 


For Case (i), we have 


dft < sup 
t>o 




dt + aidwt. 


Using the same arguments as in Theorem 14.31 yields 


Ino;* 
lim sup- 

t—^OC t 


< sup ^ ai{t) — 

t>0 


aid) 


< 0 , 


from which follows that limt_>oo Xt = 0 a.s. and the death rate of prey is exponential. 
Consider Case (ii). It follows from 

6 = ^0 + f 

do 

and Jensen’s inequality that 

E6 <Co- b[ [ Eexp{^4ds < Co - b 
Jo 

Therefore, E^t < Zt where Zt is the solution of the following differential equation 

Z't = ~b{ exp{Zt}, Zq = ^0- 

It is easy to see that Zt = — log[52^+exp{—^o}] —)• —oo as t —)• oo, then hmt_>oo E^t = 
—oo. Using Jensen inequality again gives Ext = 0. The proof of Case (iii) is similar 
to one of Theorem 14.11 We therefore omit it here. 


exp{E^s}(is. 


oi - Y - &iexp{^J 


ds 


ai{s)dv 
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